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RADIATIVE SIGNATURES OF RELATIVISTIC SHOCKS 
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ABSTRACT 

Particle-in-cell simulations of relativistic, weakly magnetized collisionless shocks show that particles 
can gain energy by repeatedly crossing the shock front. This requires scattering off self-generated 
small length-scale magnetic fluctuations. The radiative signature of this first-order Fermi acceleration 
mechanism is important for models of both the prompt and afterglow emission in gamma-ray bursts 
and depends on the strength parameter a = Xe\SB\/mc 2 of the fluctuations (A is the length-scale 
and \SB\ the magnitude of the fluctuations.) For electrons (and positrons), acceleration saturates 
when the radiative losses produced by the scattering cannot be compensated by the energy gained 
on crossing the shock. We show that this sets an upper limit on both the electron Lorentz factor: 

7 < I0 6 in /I cm~ 3 ) 7 1 / 6 and on the energy of the photons radiated during the scattering process: 

fiw max < 40Max(a, 1) (n/1 cm" 3 ) 1 ' 6 7- 1 / 6 eV, where n is the number density of the plasma and 7 
the thermal Lorentz factor of the downstream plasma, provided a < a CT n ~ 10 6 . This rules out 
'jitter' radiation on self-excited fluctuations with a < I as a source of gamma-rays, although high- 
energy photons might still be produced when the jitter photons are upscattered in an analog of the 
synchrotron self-Compton process. In fluctuations with a > 1, radiation is generated by the standard 
synchrotron mechanism, and the maximum photon energy rises linearly with a, until saturating at 
70 MeV, when a — a crit . 

Subject headings: radiation mechanisms: non-thermal — acceleration of particles — gamma rays: 
bursts 



1. INTRODUCTION 

In astrophysics, the most widely discussed mechanism 
of particle acceleration is the first-order Fermi process 
operating at collisionless shocks. It is based on the 
idea that particles undergo stochastic elastic scatterings 
both upstream and downstream of the shock front. This 
causes particles to wander across the shock repeatedly. 
On each crossing, they receive an energy boost as a re- 
sult of the relative motion of the upstream and down- 
stream plasmas. At non-relativistic shocks, scattering 
causes particles to diffuse in space, and the mechanism, 
termed 'diffusive shock acceleration', is widely thought 
to be responsible for the acceleration of cosmic rays in 
supernova remnants. At relativistic shocks, the trans- 
port process is not spatial diffusion, but the first-order 
Fermi mechanism operates nev ertheless (for reviews see 
IKirk fc Duffvlfl999 : Hillas 2005). In fact, the first ab ini- 
tio demonstrations of this process using particle-in-cell 
(PIC) simulatio ns have recently been presented for the 
relativistic case (lSt)itkovskvl l2008bl: iMartins et all [2009: 
ISironi fe Spitkovskvi^OOl l 

Several factors, such as the lifetime of the shock front, 
or its spatial extent, can limit the energy to which par- 
ticles can be accelerated in this process. However, even 
in the absence of these, acceleration will ultimately cease 
when the radiative energy losses that are inevitably asso- 
ciated with the scattering process overwhelm the energy 
gains obtained upon crossing the shock. Exactly when 
this happens depends on the details of the scattering pro- 
cess. 

In the non-relativistic case, the diffusion coefficient 
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is frequently parameterized in terms of its lower limit, 
known as Bohm diffusion, at which the particle mean free 
path equals the gyro radius. In this case, particles that 
achieve the highest possible energy radiate synchrotron 
photons up to an energy ~ 150?y (v s /c) 2 MeV, where v s 
is the speed of the shock front and n (< 1) is the inverse 
ratio of the diffusion coefficient to its value in the Bohm 
limit. Interestingly, this photon energy is independent of 
the magnetic field strength. 

In PIC simulations of weakly magnetized shocks, how- 
ever, particles appear to undergo small-angle scatterings 
on fluctuations in the electromagnetic fields that are 
driven by the Weibel instability. These have a length 
scale roughly equal to the inertial length c/uj p , where 

uj p = \J Anne 2 / '^m is the relativistic electron plasma fre- 
quency, and 7 is the Lorentz factor of the average thermal 
motion. As we show below (see Eq. ([5])) the resulting 
mean free path of a particle of Lorentz factor 7 is propor- 
tional to 7 2 , rather than the linear dependence expected 
in Bohm diffusion. This m odifies the max imum energy 
of the synchrotron photons ( Dcrishcv 2007) , but the ra- 
diation emitted when a particle is deflected through very 
small angles differs significantly from synchrotron radi- 
ation (|Landau fc Lifshit3 If 975T 1, and can, in principle, 
produce more energetic photons. This has led to the 
suggestion that "jitter" radiation from shock-accelerated 
particles is respo nsible for both the prompt emission 
ilMedvedevI 120061 ) and the afterglow (|Medvedev et alJ 
120071 ) from gamma-ray bursts. In this Letter we show 
that the inherent weakness of the scattering produced 
by Weibel-driven turbulence implies that radiation losses 
quench first-order Fermi acceleration relatively quickly. 
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In this case, the maximum photon energy is given by ([T| 
and particles are unable to produce X-ray or gamma-ray 
photons in the downstream rest frame. 

2. RELATIVISTIC SHOCK FRONTS 

Recent particle-in-cell simulations of relativistic shocks 
have revealed fundamental differences between the mag- 
netized and unmagnetized cases. The most detailed re- 
sults are available for electron-positron pair plasmas, 
which are likely to be found in gamma-ray bursts and 
pulsar winds. The upstream and downstream plasmas 
can be characterized by magnetic field strengths B u <j 
and number densities n Uj d measured in the respective rest 
frames. The upstream plasma is cold, and streams into 
the shock with a Lorentz factor 7, roughly equal to the 
thermal Lorentz factor of the downstream plasma. The 
magnetization parameters in each region can be written: 



o u = B 2 / (47rn u mc 2 ) 
<7 d = Bl/(AiT^n d mc 2 ) 



(1) 



At a parallel shock, B u = B d and <7d <C er u . If the mag- 
netic field is compressed at an oblique shock front, one 
expects er u ~ a d , since B d /B u ~ n d /n u ~ 37. However, 
if magnetic field is generat ed, then a d 3> c u . It is found 
(|Sironi fe Spitkovskyj 120091 ) that 'unmagnetized' shocks, 
with <7 U < 10~ 3 , are mediated by the Weibel instability. 
This is true also for shocks with an upstream magnetic 
field parallel to the shock normal, provided <r u < 0.1. 
On the other hand, more strongly magnetized or oblique 
shocks are mediated by the synchrotron maser instabil- 
ity that operates when particles streaming into the shock 
start to gyrate about the compressed downstream mag- 
netic field. In each case, it is possible to estimate a length 
scale on which turbulence is initially generated in the 
electromagnetic field. For Weibel-mediated (unmagne- 
tized) shocks, this is 

A w =f w c/w p (2) 

with £ w ~ 10, according to iSironi fc Spitkovskyj ([2009). 
Here, w p is the local plasma frequency. Because n<i/n u ~ 
7, the plasma frequencies in the upstream and down- 
stream media are roughly equal, and so the wavelength 
of the turbulence generated is approximately the same. 
For the purpose of estimating radiative signatures of ac- 
celerated particles, it is more convenient to characterize 
the fluctuations in terms of their 'strength', or 'wiggler' 
parameter a, defined as the ratio of their length scale 
to the length defined by the turbulent field strength: 
a = Xe\5B\/mc 2 . In the upstream and downstream me- 
dia (subscripts "u" and "d" ) one has 



a w ,d 



£ <j 1/2 

u 

6 - V2 



(3) 



This suggests that for Weibel-mediated shocks in pair 
plasmas the strength parameters are likely to be small 
(a < 1), although moderate values (a > 1) are possible 
if the self-generated field grows locally to levels close to 
equipartition. Carrying out the same analysis for un- 
magnetized electron-ion shocks results in an additional 
factor of (mi/m e ) in the value of the strength parame- 
ter, since the d ominant length- sc ale in this case is the 
ion skin-depth (|Spitkovsky||2008al ). Provided a Ui d <C 1, 
moderate strength parameters can still be expected. 



For magnetized shocks mediated by the synchrotron 
maser instability (Lyubarsky 2006), the characteristic 
length in the downstream plasma A s .d is dictated by the 
requirement that the incoming particles be significantly 
deflected: 



A s ,d =£ s ^mc 2 /eB d 



(4) 



with £ s ~ 1. This generates electromagnetic waves of 
strength parameter 



a s ,d : 



= 47 



(5) 



These waves propagate into the upstream medium with- 
out change of their (Lorentz invariant) strength param- 
eter: a su = a Sid . 

To summarize, the strength parameters associated 
with shock-generated turbulence are likely to be small 
(a < 1) in the case of Weibel-mediated pair shocks, but 
moderate a ~ 7 in the case of electron-ion shocks and 
magnetized shocks mediated by the synchrotron insta- 
bility. 

3. PARTICLE ACCELERATION 

Fermi acceleration at relativistic shocks differs from 
that at nonrelativistic shocks in that it is sensitive 
to the physic s of t he particle scattering process (e.g., 
iKirk fc Duflvl 119991 ). Typically, it is assumed that 
the particles undergo stochastic small-angle deflections 
through interactions with magnetic irregularities. Ac- 
cording to the strength parameter of these irregularities, 
the transport of the particles can then be divided into 
two distinct regimes, which we call 'ballistic' and 'heli- 
cal'. In ballistic transport, the scattering mean free path 
is shorter than the gyroradius in the local field, so that 
the curvature radius of a trajectory is larger than the 
length scale on which this quantity fluctuates. This kind 
of transport is produced by fluctuations with a < 7. On 
the other hand, in the helical transport regime, gyro mo- 
tion is only slightly perturbed. Particles have sufficient 
time to gyrate about the field while their pitch angles 
and guiding-center positions diffuse. This requires fluc- 
tuations with a > 7. This has a strong influence on 
the acceleration process and its interplay with radiation 
losses. 

3.1. Ballistic transport regime 

At a relativistic shock, a particle remains in the up- 
stream medium until it has been deflected, on average, 
through an angle of 1/7 in the upstream rest frame. A 
turbulent fluctuation of strength parameter a, deflects a 
particle of Lorentz factor 7 through an angle a/7. Pro- 
vided this angle is small, the diffusion coefficient is sim- 
ply T>e = a 2 i/ sc /j 2 , where v sc is the mean scattering fre- 
quency. The average number of scatterings in the up- 
stream medium between shock encounters is therefore 



As, 



(G) 



At each scattering, the power radiated in photons by 
the energetic particle can be estimated from Larmor's 
formula: 



A7 

7 



loss.u 



2a^ 2 7 
3mc 2 A u 



(7) 



assuming inverse Compton losses can be neglected. For 
kinematic reasons, the average energy gain per cycle is 
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roughly a factor of two (|Achterberg et al.ll2001l ). so that 
the acceleration process will saturate when the energy 
lost in the upstream medium is roughly jmc 2 . This im- 
plies iV scatt , u A7/7lio S s.u < 1, or 



7< 



/ 3mc 2 A u 7 z 



1/3 



(8) 



Applying the same argument, energy loss by scatter- 
ing in the downstream medium, where particles must be 
scattered through an angle of roughly tt/2, imposes the 
condition 



7< 



3mc 2 Ad 
2e 2 



1/3 



(9) 



which, since Ad ~ A u , is more restrictive. 

The ballistic regime requires et u .d < 7, which, accord- 
ing to (O, is fulfilled for unmagnetized pair shocks for 
all particles with 7 > 7, provided £ w < 10, and er Ui d < 1. 
For electron-ion shocks, this condition is unlikely to be 
satisfied for electrons with 7^7, but is easily satisfied 
for electrons that achieve energies comparable to those 
of thermal ions or higher, provided <7 Uj d <C 1. 

3.2. Helical transport regime 

In the helical transport regime, energy losses are im- 
portant at all points along a trajectory, and, if Bohm 
diffusion operates, the time taken to return to the shock 
front is approximately the gyro period. Under these con- 
dit ions, the maximum Lore ntz factor has been calculated 
bv lAchterberg et all (|2001l ): 



7< 



3m 2 c 3 
2e 3 B 



(10) 



At magnetized, relativistic shocks, particle accelera- 
tion by the first-order Fermi mechanism is less plau- 
sible, since, at least for superlumi nal shocks, it relies 
on st rong cross-field diffusion fe.g. JBaring fe Summerliiil 
200Sj). However, if the process does operate, particles can 
move out of the helical regime into the ballistic regime, 
as their Lorentz factor increases. Because of this, it is 
convenient to define a critical strength parameter a cr it 
such that when a = a cr j t the maximum Lorentz factor 
7max permitted by radiation losses is achieved just at 
the point at which the transport changes character from 
helical to ballistic. If a > a cr i t , all particles remain in 
the the helical regime. On the other hand, if a < a CIlt , 
particles of the maximum Lorentz factor undergo ballis- 
tic transport, but lower energy particles may be in the 
helical regime. Since the transition from helical to bal- 
listic occurs at 7 = a, the critical strength parameter in 
the downstream region follows from ((9|): 



o 2\ \ V 3 

3rac Afi \ 



2e 2 



(11) 



Expressed in these terms, the maximum Lorentz factor 

in the helical transport regime (|10p is 7 max = V a crit/ a - 
Typically, a cr j t 3> 1 in both unmagnetized and magne- 
tized shocks: inserting, for example, the length scale ap- 
propriate for Weibel-mediated shocks ([2]) one finds 



acrit = 1.4 x 10 6 4 /3 7 1/6 (n/l« 
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In the helical regime, the first-order Fermi process has 
been observed, so far, only in PIC simulation s of sublu- 
minal pair shocks (jSironi fc Spitkovskvl l2009L The fail- 
ure of the mechanism at superluminal shocks may be 
because the helical character of orbits prevents particles 
from recrossing the shock front by sweeping them away 
into the downstream medium (j Bcgclman & Kirk 119901 ; 
iBednarz fc Ostrowskil Il996l ; iLemoine et al.ll2006f ). How- 
ever, first-order Fermi is not the only possible acceler- 
ation mechanism at superluminal shocks. Strong elec- 
tromagnetic wave fields are generated at these shocks, 
and, at least in plasmas that contain some protons, these 
are thought to be responsible for part icle acceleration 
(|Amato fc Arond 120061 ; iLvubarskvi 120061 ). This can only 
happen if the nonthermal particles are confined within 
the source, i.e., if they are deflected through an angle 
~ 1 before losing their energy to radiation. In this case, 
demanding that a particle be able to complete a gyra- 
tion before radiating its energy leads to the the same 
constraint as obtained for first order Fermi acceleration 

GDp\ 

Combining the constraints from the ballistic regime © 
and the helical regime ifTU]) gives: 



for a < a cr j t 



a C rit V a crit/afor a > 



4. RADIATIVE SIGNATURES 



(13) 



(12) 



The spectrum of photons that are emitted when a par- 
ticle is scattered by turbu lent fluctuations is, in the gen - 
eral case, quite complex ([Toptvgin fc Fleishman! Il987l ) . 
However, to estimate the radiative signature of a particle 
accelerated at a relativistic shock front several simplifi- 
cations can be made. 

Consider the radiation emitted when a particle is scat- 
tered by a single fluctuation of length A and strength 
parameter a, and is unperturbed before and after in- 
teraction. The character of the emission depends cru- 
cially on the "fo rmation" or "coherence" l ength of the 
radiated photons lAkhiezer fc Shul'gal (|1987l ). At low fre- 
quencies, this length becomes large, but we confine our 
estimates to the highest frequency photons produced in 
the interaction. If a > 1, the formation length is roughly 
Acoh ~ mc 2 /eB. This is much smaller than the wave- 
length of the turbulence, so that the individual photons 
are created in regions in which the field is almost constant 
and homogeneous. The result is synchrotron radiation, 
with the emissivity defined by the local value of the field. 
The emission extends up to the roll-over frequency of the 
highest energy electrons: 

Wmax ~ Q.bj^eB/ mc 

= 0.5a 7 2 c/A fora>l (14) 

These photons are radiated into a forward directed cone 
of opening angle 1/7. Well below this frequency, the 
emissivity is proportional to cj 1 / 3 , provided the coher- 
ence length remains shorter than the wavelength of the 
turbulence. 

If, on the other hand, a < 1, the particle is deflected 
through an angle that is small compared to I/7. In this 
case, the coherence length is no longer limited by de- 
flection, but is given by the distance moved by the par- 
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Fig. 1. — The maximum Lorentz factor 7 ma x of electrons and 
the maximum photon energy £ max = 'tWmax/'nc 2 they radiate 
when scattered by magnetic fluctuations of strength a (see Eqs. J3]l 
and (O) at a rela tivis tic shock. The critical strength parameter 
a cr i t is defined in 1121 . ctf is the fine-structure constant. The jit- 
ter/synchrotron regimes are separated by the vertical a = 1 line; 
the ballistic/helical transport regimes by the a = a cr it line. See the 
electronic edition of the Journal for a color version of this figure. 

ticle in the lab. frame during the time it takes for the 
photon to move one wavelength ahead of the particle: 
A C oh ~ 7 2 c/w. The spectrum of radiated photons rolls 
over when A co h ~ A, and is flat (ocw°) towards lower fre- 
quency, as in the case of bremsstrahlung. The maximum 
frequency is roughly 



w max «0.57 I 2 nax c/A 



for a < 1 



(15) 



In each case, most of the power radiated by an indi- 
vidual electron emerges within a decade of the roll-over 
frequency that corresponds to its Lorentz factor. There- 
fore, to estimate the spectrum radiated by a power-law 
distribution of electrons, with differential number den- 
sity dn/d7 oc 7~ p , one can simply make a one-to-one 
correspondence between radiated frequency and Lorentz 
factor: u = Max(a, l)7 2 c/A. Larmor's formula, which 
holds for all values of a, states that the radiated power 
is proportional to the Lorentz factor squared, so that 
the power radiated by dn electrons of Lorentz factor 7 
is dL oc 7 2 dn. A simple change of variables then shows 
that in both the synchrotron and jitter cases, the spec- 
trum at frequencies below the roll-over frequency of the 
highest energy electrons is dL/dw oc w~ (j '~ 1 )/ 2 . 

Combining the limit on the Lorentz factor (|9]) with the 
expressions for the roll-over frequency (fT4"f and (fl"5|) , one 
finds for the maximum frequency that can be radiated 
by particles accelerated at a relativistic shock front: 



flLO-n 



( af a C rit) a < 1 

a (afa cl . it ) _1 1 < a < a c 



(16) 



where c*f = e 2 /Hc is the fine structure constant. There 
are two transport regimes: that in which the unper- 
turbed motion of the highest energy particles is ballistic 
(a < 7max) and that in which it is helical (a > 7 max ) and 



two radiation regimes: those of jitter (a < 1) and syn- 
chrotron (a > 1) radiation. This is illustrated in Fig. [TJ 

5. DISCUSSION 

Radiation emitted by relativistic electrons scattering 
in the small-scale turbulent magnetic fields generated at 
Weibel-mediated relativistic shocks has been proposed as 
the mechanism responsible for both th e prompt and af - 
terglow emission of ga mma-ray bursts (jMedvedev! 120061: 
iMedvedev et alJ 120071 1 . The evidence in favour of this 
suggestion is based on modelling the observed spec- 
tra assuming an electron distribution of power-law type 
with arbitrary high and low energy cut-offs. Power- 
law distributions are expected on theo retical grounds 
(|Achterberg et al.ll2001l : lKirk et al.ll2000h . and are indeed 
observed i n simulations of weakly magnetized, relativis- 
tic shocks (ISpitkovskvi r2008bl: ISironi fc Spitkovskvi IIM 
iMartins et al.l 120091 ) . However, the constraint on the 
maximum photon energy imposed by the above analy- 
sis (TToT) suggests that this picture is not self-consistent, 
because the scatterings are too weak to accelerate elec- 
trons to the required Lorentz factor. In order to radiate 
photons of energy ~ mc 2 in the plasma rest frame, strong 
fluctuations of large length-scale with a ~ afa crit ~ 10 
are required. 

The above conclusion rests on the assumption that 
the same fluctuations are responsible for both the parti- 
cle transport and radiation. In terms of the strength 
parameter a and length scale A that we use to char- 
acterize the fluctuations, the deflection angle scales as 
Ad oc a and the radiation losses as A7 oc a 2 /A. If, 
therefore, the scattering responsible for isotropization oc- 
curs on fluctuations of comparable strength, but much 
larger length scale than those responsible for the radi- 
ation losses, the limit on the maximum photon energy 
(|16|) is relaxed. In principle, the fluctuations induced 
by the Weibel instability could be responsible for pho- 
ton production, provided longer wavelength fluctuations 
are present to provide the necessary isotropization and 
transport. The accelerated particles themselves appear 
to generate l onger wavelength fl uctuations downstream 
of the shock (|Keshet et al.l 120091 ) . but this is a relatively 
small effect compared to that needed to significantly in- 
fluence the maximum photon energy. 

On the other hand, if, as simulations suggest, the 
Weibel-induced fluctuations are responsible for the trans- 
port, the bulk of the radiation must be produced by in- 
teraction with fluctuations of much shorter wavelength. 
An obvious candidate is the soft photon field pro- 
duced by the interaction of thermal electrons with the 
Weibel-induced fluctuations — the 'jitter' analog of the 
synchrotr on photons produced by relativistic thermal 
electr ons (Baring & Braby 2004: iGiannios fc Spitkovskvi 
.2009). With these photons as targets, the radia- 
tion mechanism is analogous to the synchrotron self- 
Compton mechanism, which has been discussed in con- 
nection with the problem of rapidly decaying mag- 
netic fluctuations (jRossi fc Reesl 1mm . While this 
model offers a simple explanation for the produc- 
tion of high energy photons, it does not easily ac- 
commodate those low- frequency BATSE burst spectra 
(|Baring fc Brabvi 120041 ) that violate the s o-called syn- 
chrotron 'line of death' ([Preece eTaT] [19981 ). These rare 
cases might be accounted for if the upscattering occurs 
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deep in the Klein-N ishina regime (Deri shev et al.l [20011 ; 
iBosniak et al.l [2009h . an explanation which is facilitated 
if the target photons arise not from synchrotron emission, 
but from relatively hard jitter radiation. 
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